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1. Introduction 

The measurement of the temperature-dependent magnetic susceptibility, x(T), pro- 
vides a standard essential diagnostic method for establishing the magnetic properties of 
a system. A careful comparison between measured data and the predictions for x(T) 
as derived from a model Hamiltonian is routinely performed with the goal of estab- 
lishing numerical values of model parameters, for example, the exchange constant (s) of 
the Heisenberg model of interacting isolated spins. The success of this technique is so 
firmly established that it is taken for granted that there is a one-to-one correspondence 
between a given form of x{T) and the numerical values of the model parameters. Most 
certainly it is unnatural to contemplate that one might be able to continuously vary 
the parameters of a model Hamiltonian and yet generate a single, invariant form for 
x{T) and similarly for other thermodynamic quantities. Yet surprisingly, there are a 
number of exceptional systems, where there is a continuous infinity-to-one correspon- 
dence between model Hamiltonians and measurable thermodynamic quantities. One of 
these exceptional cases was recently encountered in the course of attempting to 

determine the exchange constants of a simulating Heisenberg model from experimental 
susceptibility data for a specific synthetic magnetic compound |J. 

We shall refer to systems of a continuous family having the same eigenvalue spec- 
trum as being "isospectral" . This notion is chosen in analogy to the use of "isospectral- 
ity" in other areas of physics, e. g. the occurrence of supersymmetric pairs of Hamilto- 
nians (see, e. g. [Q, 7.1) or the problem of bounded domains with isospectral Laplacians 
("Can one hear the shape of a drum ?" ||). The subject of isospectral spin systems is 
not completely novel but has been discussed in the literature only on a few occasions, 
e - g- @) [0- However, there is, to our best knowledge, no systematic account of this 
phenomenon, the first steps of which will be presented in this paper. 

In addition to providing a general approach to isospectrality we analyze in depth 
two cases of distinct Heisenberg systems where continuous variation of the exchange 
constants gives rise to one and the same set of temperature-dependent thermodynamic 
quantities. The operational conclusion for an experimentalist is quite sobering, in that 
for these specific systems measured data alone cannot fix the exchange constants. Com- 
parison between theory and experiment can only place a weak constraint on a continuous 
family of equally acceptable choices of parameters. Although we provide some helpful 
insights, it is very difficult to formulate the general set of conditions to be met so as 
to achieve such exceptional model systems. It is reasonable to expect, that if one and 
the same temperature dependent thermodynamic quantity is generated by a continuous 
family of Hamiltonians then necessarily all members of that family share the very same 
eigenvalue spectrum. This is indeed the case as it is proven in section 3 for two partic- 
ular thermodynamic functions. 
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In short, our goal in the present work is to provide a first systematic study of 
Heisenberg isospectral spin systems. The remarkable advances 0,|§ in synthesis mag- 
netochemistry, of incorporating significant numbers of interacting paramagnetic centers 
within individual molecules, may provide the impetus for wider studies that will yield 
a more comprehensive set of conditions for the occurence of isospectral spin systems. 

The layout of this paper is as follows: 
In section 2 we introduce our notation and the basic concepts of "isospectrality" , 
"complete isospectrality", and "covariant isospectrality" for spin systems with 
Heisenberg Hamiltonians. Families of isospectral systems are algebraic varieties in 
the space J of coupling constants. Covariant isospectrality is implemented by a 
unitary representation of some Lie group, which simplifies the calculations considerably. 
Unfortunately this is a rare case, as we will see. Complete isospectrality means that 
all eigenvalues of two systems with the same magnetic quantum number are in 1:1 
correspondence and equal. We do not know whether this is a strictly stronger property 
than plain isospectrality, except for the case of a trivial counter-example. However, 
we need this apparently stronger concept to derive the conclusion that completely 
isospectral systems share the same magnetic susceptibility function. This is done in 
section 3 where also the inverse problem is settled as well as the analogous question for 
the specific heat function. The result in short is the following: Plain isospectrality is 
equivalent to possessing the same specific heat function and necessary for possessing 
the same magnetic susceptibility function. Complete isospectrality is sufficient for 
possessing the same magnetic susceptibility function. In section 4 we identify the 
isospectral invariants which are linear or quadratic in the coupling constants. This is 
crucial for section 5 where we show that the triangle with spin s — | is both of completely 
and covariantly isospectral type but that for s > | isospectrality breaks down. On the 
other side, if the number N of spin sites exceeds 3, covariant isospectrality is no longer 
possible. This is proved in section 6 with the aid of MATHEMATICA® 4.0 and some 
trace formulae which are explained in Appendix A. The tetrahedron (N = 4) with s = \ 
nevertheless possesses completely isospectral families of dimension one and two, as is 
shown in section 7 and Appendix B. In section 8 we provide some heuristic arguments 
in order to explain the findings of the previous sections. We expect that isospectrality 
only occurs if the number of all possible bonds exceeds the number of independent 
eigenvalues and show that this never happens except for s = | and N = 3,4,5. A table 
summarizing our results and conjectures on the occurence of isospectrality for different 
N and s and concluding remarks are provided in section 9. 

2. Notations and definitions 

We consider spin systems with N spin sites, spin quantum number s and isotropic 
Heisenberg coupling between all sites x and y with coupling constants J xy . For sake of 
compact notation we will write the (^) coupling constants J xy as the components of a 
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vector J E J = IRv 2 J. Thus a specific point of J uniquely specifies the strength of the 
interactions between all pairs of spins and will be sometimes called a "system" . 

Let Sx \ (i = 1, 2, 3), denote the three components of the spin observable S x at site 
x and, as usual, 

s = s x , s« = s ± = sw ± is®, (i) 

X X 

denote the total spin vector and its various components. All linear operators occuring 
in this context will be identified with the corresponding dim x dim-matrices, dim = 
(2s + 1)^ being the dimension of the total Hilbert space of the spin system, w. r. t. the 
fixed basis consisting of tensor products of eigenvectors of ■ The Hamilton operator 
can then be written as 

H = J ■ H = ^2 JxyHxy, (2) 

x<y 

where 

3 

H xy = S x -S y = Y,S<!>S<p. (3) 
i=\ 

Here H is an (^) -dimensional vector the components of which are dim x dim- 
dimensional matrices H xy . 

If the spin system is coupled to a constant external magnetic field H, the total 
Hamilton operator will be 

H(h) = H -hS {3) , (4) 

where h = g^B^ contains the common combination of the gyromagnetic ratio g and 
the Bohr magneton /ib- As usual, the partition function, which yields all the standard 
thermodynamic functions, is defined by 

Z(p,h)=Tr(eM-pH(h)))- (5) 
In particular, one obtains from Z the specific heat function 

c((3)=(3 2 -^\nZ((3,0), (6) 

and the magnetic zero field susceptibility function 

X (f3) = ^^Z(P,h)\ h=0 . (7) 

Further we will need the traces of powers of H , t n = Tr(i?o ), n = 0, 1, 2, ... , dim, and 
the set of traces 

T(H ) = {t n \n = 0,1,2,..., dim}. (8) 
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Two Hamilton operators with the same N and s, = J (1) • H and = ■ H 
are called isospectral if they have the same eigenvalues, counted with multiplicity, or, 
equivalently, if they generate the same characteristic polynomial: 

det{H^ 1] - A) = det(^ 2) - A) VA e R. (9) 

According to the above remarks, we will also speak of "isospectral systems" . Sometimes 
we will apply the term "isospectral" to more general pairs of operators derived from 
H®,i = 1,2 if there is no danger of confusion. Clearly, (P) defines an equivalence rela- 
tion ~ on J . The coefficients of a characteristic polynomial P(X) = det(H — A) can be 
viewed as polynomials of the coupling constants J xy . These polynomials assume constant 
values exactly on the ~-equivalence classes [J]~, which we will call isospectral classes. 
Consequently, the isospectral classes are algebraic varieties in J', since they are defined 
by a finite number k of polynomial equations p v {J) — 0, v = 1, . . . , k^ 



In algebraic geometry there are various equivalent definitions of the dimension of 
an algebraic variety, which are, however, too technical to be reproduced here (see, for 
example, [|l(|, chapter 9). For our purposes it will suffice to note that for the special case 
where the Jacobian matrix ( w- 1 has locally a constant rank r, the corre- 

V 3 J v=l...k,j=l...( „ | 



-r. 



sponding isospectral classes are locally differentiable manifolds of dimension £ — ( ) 
This is an immediate consequence of the fibration theorem (see, for example, JI2(] , the- 
orem 3.5.18). Especially, if k > ( 2 ) and the rank of the Jacobian is maximal, r = 
the isospectrality classes will only consist of a discrete point set. Note, however, that 
the equivalence classes will never be trivial since they at least consist of the orbits of 
the group of discrete symmetries of J generated by permutations of spin sites. 



We will say that the pair (N, s) is of isospectral type if the corresponding space 
of coupling constants J contains at least one isospectral equivalence class of dimension 
i > 1. The largest dimension I of isospectral equivalence classes will be called the 
dimension of the isospectral type. 

Functions / : J — > K which are constant on isospectral classes will be called 
isospectral invariants. 



Moreover, we will consider a special case of isospectrality in which the equivalence 
classes are easily calculated. Obviously, two isospectral Hamiltonians Hq and i^ 2 '' can 
be related by a unitary transformation U : 

Hf ] = U*H^U. (10) 

j For an elementary introduction into the theory of algebraic varieties, see e. g. |fi0[ . Note that 
algebraic varieties need not be differentiable manifolds since they may contain "boundaries" like the 
vertex of the light cone x 2 + y 2 + z 2 = c 2 t 2 or "hairs" like x = y = in the variety (x 2 + y 2 )z = 0. 
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U maps the eigenvectors of onto the corresponding eigenvectors of Hq. In the 
case of a one-dimensional isospectral equivalence class parametrized by a coordinate t, 
U can be chosen to depend smoothly on t. Following a closed loop, the corresponding 
unitary transformation need not reduce to the identity transformation, but may include 
some phases. We will come back to this phenomenon later. If for sufficiently many 
curves of isospectrally equivalent points the corresponding unitary transformations U(t) 
can, moreover, be chosen to be a one-parameter group, the system (N, s) will be called 
of "covariant isospectral" type. More precisely, we define: 

Definition 1 (N, s) is of covariant isospectral type iff it is of isospectral type and for 
any isospectrally equivalent which are sufficiently close, there exists an anti- 

Hermitean dimx dim-matrix fi, and a real ( 2 ) x ( 2 ) -matrix M and some to G R, such 
that for allteR 

exp(n*t)J {1} -Hexp(nt) = (exp(-Mt) J (1) ) • H (11) 

and 

exp(-Mt )J (1) = J (2) . (12) 

It follows immediately, that the orbit exp(— Mi) in ^7-space corresponds to a family 
of isospectral Hamilton matrices. Since in our case all Hamilton matrices are real and 
symmetric, Q can also be chosen as real and hence anti-symmetric, and thus exp(fit) 
represents a rotation in the real Hilbert space M. dvm . 

The condition (|TT| ) may be replaced by its equivalent infinitesimal version: 

[H,n] = -M*H. (13) 

This equation could be used to show that M, considered as a linear transformation in 
the space of (^) -dimensional vectors with matrix entries endowed with the scalar prod- 
uct Ti(K ■ H), will be an antisymmetric matrix. We will give an independent proof of 
this fact later. 



Further we note that the set of solutions (Q, M) of ( p~3|) will be a Lie algebra with 
respect to the obvious vector and commutator operations. Hence covariant isospectral 
equivalence classes will be orbits of the corresponding matrix Lie groups. 

Thus far we have only discussed isospectral systems in the absence of an external 
magnetic field. If we include the field TC and allow for a corresponding /i-dependence 
of the Hamiltonian, we have to consider a slightly stronger concept of "complete 
isospectrality" . 

Definition 2 Two Hamiltonians Hq and (or, equivalently, two vectors and 
j(2) £ j) are ca n e( i completely isospectral (J^ ss J^) iff for all h G R H^(h) = 
- hS<® and H (2 \h) = - hS® are isospectral. 
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Proposition 1 The following conditions are equivalent: 

(i) J«^J( 2 ) ; 

(ii) P M H^P M ~ Pm^ 2) Pm 

for all projectors P m onto the eigenspaces of 
corresponding to the eigenvalue M , 
(Hi) FsH^Fs-FsH^Fs 

for all projectors F$ onto the eigenspaces of S 2 
corresponding to the eigenvalue S(S + 1). 

Proof: (i)=>- (ii). Since all H(h),h € R, commute, there exists a system of joint 
eigenprojectors P^m,, such that 

H{h) = Y,^v-hM v )V vMv . (14) 



Hence 

P MJ Ff(/i)P M = ^ {e v -hM v )V vMv . (15) 

i/,ATi, ,M„=M 

If ~ H^\h) for all heR, then both systems have the same set of eigenvalues 

e v — hM v with the same multiplicities TrV v ,M v - Hence, by (jlBD , also Pjyii'( 1 '(/i)P A f ~ 
P m #( 2 )(/i)Pm, especially P M ^ 1} P M ~ P M H®P M for all M. 

(ii) =^(iii). The eigenspaces corresponding to V Vi m v can be further split into 
eigenspaces of S 2 . This can be written as 

Smax 

V v ,m v = E ^v,M v ,s- (16) 

S=\M U \ 

Applying the ladder operators S gives 

Tr V V)MvtS = Tr V v , s , s = Tr V v , s for all M„ = —S, ...,S. (17) 

These numbers are the multiplicities of the eigenvalues of the operators FsP mH^FsP m, 
which are hence isospectral for i — 1,2. By summation over M we conclude that also 
FsH^Fs, are isospectral for % = 1,2. 

(iii) =^(ii). This can be shown analogously by considering 

F S H F S = e » V U> ( 18 ) 

v,S„,S„=S 

s v 

Ks v = E V U,M (19) 



M=-S v 



and 



TrP^ = (2^ + l)Tr^ M . (20) 
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(ii)=^(i). This follows from 

H{h) = ^{P M H P M -hM). 



(21) 



M 



Obviously, ~ implies ~ J( 2 \ but not conversely: 

/ 1 \ / -1 \ 



Take N — 3,s — I, then J« 



1 

V 1 / 



J(2) 



, but 96 J 



(2) 



v- 1 / 



Unfortunately we do not know of less trivial counter-examples. The problem is the 
following: 

Being rotationally symmetric, Hq commutes with S 2 and S^ 3 \ hence each eigenspace 
of H with eigenvalue e v is spanned by simultaneous eigenvectors of S 2 and say, 
\u, X,/jl), v — 0, . . . , D, A = 1, . . . , d v ,fi — —S„ t \, . . . , S Vj \, such that 

5 2 |z/,A,/i) = S v , x (S v , x + l)\iy,X,fj), 5 (3 V,A,/i) = /i|z/,A,/i). (22) 

The degeneracy of the eigenvalues e v will be 



n v = ]T(2^ A + 1). 



(23) 



A=l 



We do not assume that the S^a have different values. The corresponding eigenvalues of 
H(h) =H - hS® are 

E v ,\,n = e v - hpi. (24) 

For "generic" J G J we expect that H will have no further degeneracy besides that 
dictated by rotational symmetry ("minimal" degeneracy), i. e. we expect that d u = 1. 
So A can be skipped and \i = —S u , . . . , S u . In this case, "isospectrality" and "complete 
isospectrality" would be equivalent. So, heuristically, we may consider these two notions 
as having equal meaning, although we have to distinguish between them for the sake of 
mathematical rigour. 

For the most important case of isospectral systems which are obtained by 
continuously varying the coupling constants we can show the following: 

Proposition 2 Any two systems joined by a continuous curve of isospectral systems 
are completely isospectral. 

Proof: Consider a curve t 1— > J(t) and consider 

H (t) = H-J(t) = J2^(t) (25) 



and 



(26) 



Since 1 1— > V u (t) is continuous and Tr(V u (t)Fs) assumes only non-negative integer values, 
the latter must be constant w. r. t. the parameter t. However, Tr(V u (t)Fs) equals the 
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multiplicity of the eigenvalue e u in (|26|). Thus all FsH (t)Fs, t G M, are isospectral, and, 
by proposition [I], all i?o(^) are completely isospectral. ■ 
Note that, according to this proposition, in the above counter-example the two 
systems cannot lie in the same connected component of an isospectral class. 



We add some definitions concerning symmetrical polynomials which will be of later 

use: 

The equation 

d d 

Y[(x + x n ) = J2^x d -" (27) 

n=l u=0 

defines the elementary symmetrical polynomials 

d d 

si = ^x n , s 2 =^2x n x m , ... s d = JJx n . (28) 

n=l n<m n=l 

These also appear, up to a sign, as the coefficients of the characteristic polynomial of 

d 

p(X) = l[(X-E l ), (29) 

i=l 

written as polynomials of the Ei, where d = dim. Every other symmetric polynomial 
can uniquely be written as a polynomial of the s v (see, for example [|IIJ], 7.1, Theorem 
3). This holds especially for 

d 

t n = Tr(HZ) = Y,E?. (30) 

i=l 

Conversely, each s u can uniquely be written as a polynomial of the t n ,n = 1 . . . d, (see, 
for example UlO |, 7.1, Theorem 8) e. g. 

s 5 = — (t\ - m\t 2 + lbhtj + 20^3 - 20^3 - 30tit 4 + 24t 5 ) . (31) 
5! 

This representation is independent of the dimension d. 



3. Specific heat and magnetic susceptibility 

As mentioned in the introduction, (completely) isospectral spin systems will give rise to 
the same thermodynamic functions like specific heat and magnetic susceptibility. In this 
section we will state this more precisely and also prove the converse, up to the subtle 
distinction between complete and plain isospectrality. 

Lemma 1 Two spin systems are isospectral iff T(H^) = T(Hq ). 

Proof: Recall that T(H ) was defined as the set of traces t n = Ty(Hq), n = 0, . . . , dim. 
Hence the "only if part is obvious. From the remarks at the end of section 2 it follows 
that if two systems possess the same t n ,n e N, they share also the same values of 
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the standard symmetric polynomials s u (E) and hence have the same characteristic 
polynomial p(X). I 

Proposition 3 Two spin systems are isospectral iff they possess the same specific heat 
function. 

Proof: c(/3) can be expanded into a Taylor series at /3 — 0: 

W = (-| + ~) f (32) 

+ .... 



This is the starting point of the so-called "moment expansion method", see e. g. |IT| 7.3. 
Obviously, each coefficient of (3 n uniquely determines t n , if the other traces t m ,m < n, 
are already known. Note that t = dim, t\ = 0. Together with lemma [1] this completes 
the proof. ■ 
Now we consider again H{h) = Hq — hS^ with eigenvectors \u, A,/i) according to 
the previous section. By its very definition, two completely isospectral systems share 
the same partition function Z{f3, h) and any other thermodynamic function which can 
be derived from it. Especially, the following holds: 

Proposition 4 Two completely isospectral systems possess the same magnetic 
susceptibility function. 

To tackle the converse problem, we consider Z(/3, 0) = J2 U n u exp(—j3e u ) and define 
the coefficients a v implicitely by 

Tr (S (3)2 exp(-/?# )) = a » n » e M~^) (33) 

V 

In the case of minimal degeneracy, i. e. d v = 1, we have n v = 2S V + 1 and 

s, 



^7 E ^ = \s»{Sv + l). (34) 

In the general case, 



d v 

1 

0~u 



A=l n=—S Vt \ 
d, 



- f2^S u>x (S u>x + l)(2S v>x + l). 



y A=l 



Proposition 5 Two spin systems with the same susceptibility function are isospectral. 



Isospectral spin systems 11 
Proof: Since 

X(P) = Tr (S^ exp(-/Hf )) , (36) 

the two systems will have the same function 

_ Tr (eM-(3Ho)S^) 
W3) ~ Tr(exp(-^ )) (3?) 
E?=o <r v n v exp(-/3e v ) 
Ef =0 ra„exp(-/3e„) 

Since e < £1 < £2 < • • •, the terms exp(— /3e„) are of different orders of magnitude for 
f3 ^ oo. The first term increasingly dominates, hence 

a n exp(-/3eo) , QS ^ 

hm f(P) = = (Jo- 38 

/3->oo n exp(-/3£ ) 

If we subtract this limit from the dominant term asymptotically becomes 

n exp(-/5e ) 

^iexp(-/3(ei-e )). 
n 

In the next step we have 

f{(3) -a -— exp(-/?( ei - e„)) —2 exp(-/3(5 2 - e )), (40) 

n n 

and so on. In this way, from the behaviour of f(/3) for (3 — > oo, we may extract the 

values 

a o, • • • • • • and £l - e , e 2 - e 0> ...,£„- e , ej3 - £o . 

Let t n = Tt(Hq) as above and u n = Tr(#™S( 3 ) 2 ), n e N. u = Tr^ 3 ) 2 ) can be 
calculated independently of H . Since 

D D 

fMi = Y]^un v = noY]- ! ^-, (41) 

i/=0 f=0 u 

no and hence o^n^, z/ = 1 . . . D, are also uniquely determined. 
Next we consider the Taylor expansion of f(/3) at (5 — 0: 

/C9) = ^_a /J+ (_^! + ^y (42) 

+ ^°^4t 3 24%) 6t 2 Atl + 2AtJ P 
+ ... 

Recall that t , /i are known. The linear term then gives 

D 

IH = Tt(H S^ 2 ) = ^n v e v . (43) 

vo=0 
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On the other side we know the 1. h. s. of 

D 

22 Ovn v {e v - e ) =111- EqUo, (44) 

i/=0 

hence Eo and all £ v , v = 1 . . . D are also known. Similarly, 

D 

22 °vn v {E v - e ) 2 =Hz- 2e jt/i + elfi , (45) 

z/=0 

hence fi^ is known and from the /3 2 -term in (|^) also t 25 and so on. 
Eventually, we obtain all t n , /i n , n = 2 . . . dim solely from xifi)- According to lemma [I] 
this gives us all eigenvalues of H with multiplicity, i. e. n u , v = 0, . . . D and the two 
spin systems are isospectral. ■ 
The proof does not give complete isospectrality: If some eigenvalue e v belongs to 
different 3 vX , A > 2 then from n v = J2x ( 2 ^a + 1) and a v = Zx ^^jffi^ the 
S U) \ cannot be uniquely determined. 

4. Some isospectral invariants 

Criteria for non-isospectrality could, in principle, be checked by brute force methods: 
Calculate the characteristic polynomial of the matrix H = JH, say p(X) = J2t=o c v^ v ■ 
Select (^) different coefficients c$, cj, ... [y = dim being excluded since c^ m = 1) and 
calculate the Jacobian determinant 

Jac(J) = d ^ Ci "y\ (46) 

preferably by using a computer algebra software. If the Jacobian nowhere vanishes, 
according to remarks in section 2, (N, s) cannot be of isospectral type. 

In practice this method will, even for small N and s, rapidly become extremely 
memory- and time-consuming. To simplify the problem one could - in the case of 
complete isospectrality - restrict oneself to subspaces invariant under H , for example 
subspaces TC(M) of constant magnetic quantum number M. 

The space with maximal M, TC(M = Ns) is one- dimensional and is spanned by the 
product state 

y?o = \s,s, ...s) (47) 
which is an eigenstate of J ■ H for all J G J with eigenvalue 

E = s 2 J = s 2 Y,J*y (48) 

x<y 

This proves the first part of the following 
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Lemma 2 // J^ 1 ' • H and ■ H are completely or covariantly isospectral, then 

£4i> = £4?, («) 

z. e. J 1 - 1 ) and «j( 2 ) Zze in the same hyperplane perpendicular to 1 = (1, 1, ... 1). 

Proof: If all J xy > 0, then £7 — s 2 J w iU be the maximal eigenvalue of ■ H. 
In fact, ((po\H xy (po) is the maximal expectation value for each H xy . Hence J is an 
isospectral invariant at least in the domain J + = {J\ all J xy > 0}. Now assume 
covariant isospectrality and let t \— > J{t) = exp(— Mt)J(O) be an isospectral curve 
which will be restricted to J + . According to what has been said before, = J = 
|(J|1) = ( J|l) = (-MJ\1) = (J\ - M*l) = (J\M1) for all J G J+. Since J+ 
linearly generates J it follows that Ml = 0, i. e. all row sums of M vanish and J is an 
isospectral invariant on the whole space J . M 
Before preceding with M = Ns — 1 we will show that also 

Lemma 3 ||J|| 2 = ^2 x<y J xy is an isospectral invariant. 

Proof: Obviously, Tt(Hq), the sum of all eigenvalues squared, is the same for isospec- 
tral Hamiltonians. After expanding the square (J ■ H) 2 one realizes that only those 
products H xy H uv have a non-zero trace where x = u and y = v. Hence Tt(Hq) = 
E x<y Jly^iHly) = (E x<y J 2 xy ) ■ + 1) 2 (2* + 1)*, see Appendix A. The actual 

value for Tr(H 2 y ) is irrelevant for the proof; what matters only is that it is independent 
of x, y. This concludes the proof. ■ 

From Lemma ||| it follows immediately that in the covariant isospectral case the 
matrix exp(— tM) leaves the Euclidean norm ||...|| invariant and thus must be an 
orthogonal transformation and its generator — M will be antisymmetric. 

5. The triangle (N = 3) 

5.1. 8 = 1 

Next we consider eigenvalues with eigenvectors in the subspace TC(M = Ns — 1), but 
restricted to the case N = 3. With the abbreviations 

J = Jl2 + J23 + Jl3, (50) 
r = J12J23 + J12J1Z + ^23^13; (51) 



the eigenvalues are calculated to be (see [13|], § 62, Exercise 2 ) 



E = s 2 J, E h2 = s(s - 1) J ± VJ 2 - ST. (52) 

Hence the first three eigenvalues are constant on curves with constant J and T. These 
are circles with radius 



^J 2 -2r, (53) 



Isospectral spin systems 



14 



the center of which is located on the line J i2 = J23 = J13, including the degenerate case 
r = 0. 



For s — ~ the list of eigenvalues is already exhausted: Due to rotational symmetry 
the value E is 4-fold degenerate (S = §) and the E 1<2 are two-fold degenerate (S = \). 
We conclude: 

Proposition 6 The system N = 3, s = ~ is of complete isospectral type with dimension 
2. 

We now consider the question whether the triangle with s — is of covariant 
isospectral type, i. e. we seek solutions of 

[Jf , Q,) = MH, (54) 

Q and M being anti-symmetric. Let T = T 3 be the unitary left shift operator which 
represents a cyclic permutation of the spin sites. Then a solution of ( |54| ) is given by 

/ 1 -1 \ 



n=-^=(T-T*), M=^-= 
2^ V3 



\ 



-1 1 
1 -1 



/ 



0, can also be written as 

i 



4V3 



cr 3 • (<ri x a 



2)- 



(55) 



(56) 



where the crj, {i = 1, 2, 3) denote the Pauli-matrices. Obviously, Q is rotationally 
symmetric which entails complete isospectrality. The factor is chosen such that the 
parameter t in exp(tM) will be just the angle of rotation. T 3 = T entails Q 3 = — 10, 
hence the exponential series of exp(tfi) will be actually a polynomial in Q: 



exp(tfi) = 1 + 2tt sin - + (2fi) 2 (l - cos ■ 
For special values of t we obtain: 

47T 

exp( —Q) = T, 
exp(47rfi) =T 3 =1 
exp(27rfi) = jj(T 2 + T-l). 
The last expression can be rewritten using 



(57) 



(58) 
(59) 
(60) 



H 



1 

V 1 / 



H = - Ps - Pi 
4 V 2 2 



where P3 (resp. Pi) denotes the projector onto the subspace S = | (resp. S = | 
result is 



exp(2vrfi) = Ps - Pi. 



(61) 
The 
(62) 



Isospectral spin systems 



15 



This means that an eigenstate of H with S = | acquires a phase of it after a full 
rotation in J^-space, analogous to the occurence of Berry phases for adiabatic loops in 
parameter space. Summarizing, we state the following proposition which, in essence, is 
due originally to V. G. Grachev 0: 

Proposition 7 The system N = 3, s = 1/2 is of both completely and covariantly 
isospectral type. 



5.2. A physical example 

An interesting and timely application of this theory is provided by the example of the 
molecular magnet (CN 3 H 6 )4Na 2 [H4V 6 P403o(CH 2 ) 3 CCH 2 OH 2 ]« 14H 2 which features 
two uncoupled systems of three V 4+ (s = |) ions that interact via antiferromagnetic 
Heisenberg exchange. It has been proposed J| that the Coulomb interaction between an 
Na ion and two of the three V A+ ions gives rise to what is essentially an isosceles triangle, 
with the distances between the three vanadium ions being 3.20, 3.21, and 3.36A. It is then 
quite reasonable to assume, that the three exchange constants satisfy J i2 = J i3 ^ J 23 • 
In fact, calculation of the weak field susceptibilty has yielded results that are in excellent 
agreement with accurate susceptibility measurements from room temperature down to 
2K upon assigning the values Ji 2 = J13 = 64. 7K and J23 = 7.5K 0] |], [0. Moreover, 
the calculated energy level spacings that follow from these assignments have recently 
been confirmed to good accuracy in a direct manner by inelastic neutron scattering [14[ . 



Nevertheless, as the work of this section has shown, the identical energy levels and the 
identical temperature dependent susceptibility emerge for the continuous choices of the 
three different exchange constants that lie on curves with J = 136. 9K and T = 5156. 6K 2 



5.3. s>\ 

For s > I we consider the next subspace 7i(M = 3s — 2). The characteristic 

polynomial p(A) = Ylt=i c ^ ^ as Deen calculated using MATHEMATICA® 4 .0, but is 
too complicated to be presented here. One particular Jacobian reads 

d(c 3 ,C 4 ,C 5 ) = ^ _ _ _ gs + 6s 2 )(1 _ g s + 15s2) _ (63) 

oyJ\2i ^23, Jvz) 

This is a polynomial in s which has no integer or half-integer roots. Therefore we have 
proved the following 

Proposition 8 The system N = 3, s > | is not of complete isospectral type. 

% It is usual to measure the coupling constants in units of Kelvin. The corresponding energies are 
obtained by multiplying with the Boltzmann constant ks- 
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6. Isospectrality for N > 3 

The question arises whether our result that complete isospectrality only occurs for s — ^ 
also holds for N > 3. Our method of calculating the Jacobian ([iffi) for arbitrary s will 
no longer work for larger N. However, we can prove a weaker statement, namely 

Proposition 9 Systems with N > 3 cannot be of covariantly isospectral type. 

Proof: For this we need a trace formula which will be explained in the Appendix A: 

Tr(# 3 ) = 5>^) 3 (-l s2 ( s + l ?( 2s + V") ( 64 ) 

x<y ^ ' 



x<y<z 

From the isospectral invariance of Tr(H 3 ) we conclude that also 

f 3 = J2(Jxy) 3 -Hs + l) Yl J *y J yz J *z ( 65 ) 

x<y x<y<z 

will be invariant. Now we consider four different spin sites (using N > 3) denoted by 
1,2,3,4 and consider vectors J G J which have vanishing components except possibly 
for J i2, J13, J23, J14,, J24, -hi- One-parameter isospectral curves passing through J satisfy 

J t h(J(t)) = 0. (66) 

Using © and = -MJ(t), © can be written as an equation which is linear 

w. r. t. the 15 relevant matrix entries of M and trilinear w. r. t. the 6 non-vanishing 
components of J. Using MATHEMATICA® 4 .0 it is easy to show that (|66|) has only 
the trivial solution M — 0. For example, one may randomly choose 15 vectors J (it 
suffices to consider components —1,0,1) and cast the corresponding 15 equations of 
the form (^) into matrix form. Non-trivial solutions exist only if the determinant of 
this matrix, which is a polynomial in s of degree 30 with integer coefficients, vanishes. 
However, the zeros of the polynomial can be numerically computed and shown not to 
attain half-integer or integer values. I 



So it seems that the concept of covariant isospectrality is of little use having only 
one single application for N = 3,s = 1/2. However, covariance may be restored for 
iV > 4 if the class of admissible Hamiltonians is suitably extended, e. g. to include also 
Hamiltonians which are bi-quadratic in the spin observables. However, this is beyond 
the scope of the present article. 

Of course, our proposition |9| does not exclude plain isospectrality for N > 3. 
Indeed, we will show that the system N = 4, s = | is completely isospectral, albeit 
not covariantly isospectral, in the next section. 
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In the case N = 4, s = 1/2 it is still possible to calculate the coefficients of the 
characteristic polynomials of H restricted to the subspaces with M = 0, S = 2,1,0. 
Obviously, this is enough in order to study complete isospectrality since all eigenvalues 
of Hq appear within these subspaces. It turns out that all coefficients can be written 
as functions of four fundamental invariants Ii, I 2 , 1 3 , 1^. These can most conveniently 
be written in terms of new coordinates in J , which are defined as half the sums and 
differences of the coupling constants of adjacent edges: 

S\i = + ^34), ^12 = ^(^12-^34), (67) 

etc. 



Proposition 10 Two spin systems with coupling constants J^,J^ with N = 4, s = 

1/2 are completely isospectral iff the following four functions assume the same values 
for J« and : 

h = D 2 2 + D 2 3 + £>? 4 , (68) 

1 2 = S 2 l2 + S 2 3 + S 2 14 , (69) 

1 3 = S 12 + S 13 + S14, (70) 

h = 2D 12 D n D u + D\ 2 S l2 + D 2 3 S 13 + D 2 U S 14 - S 12 S 13 S 14 . (71) 

Now let I' be the functional matrix obtained by partial differentiation of Ji, I 2 , 1%, I4 
with respect to its 6 arguments Si 2 , • ■ ■ , -D14. The rank of I' assumes its maximal value 
of 4 iff no determinant of the 15 possible 4x4 submatrices of /' vanishes. We denote 
the subset of those points with maximal rank by 1Z C J . If an isospectral class lies 
entirely within 1Z it will be a 2-dimensional submanifold of the 6-dimensional space J . 



This follows by a well-known theorem of differential geometry (see e. g. |T3|] Theorem 
3.5.4). We will call this case generic, the other cases exceptional. 



Although J is six-dimensional, one can visualize the 2-dimensional submanifolds 
in the generic case. 1% = const, defines a sphere in the 3-dimensional D-sp&ce with co- 
ordinates -D12, -D13, -D14. I2 = const, and I3 = const, define the intersection of a sphere 
and a plane, i. e. a circle in S-space. For given D the last equation Ia(D, S) = const. 
picks out a finite number (actually < 6) of points in that circle. If the corresponding 
angles are drawn as different radii r u = 1 + ^ in .D-space (identifying points with 
r v = 1/2 and 3/2) we obtain a surface folded in a complicated way. 



A large number of 2-dimensional generic isospectral classes have been identified 
numerically. The exceptional classes are one- or zero-dimensional and will be discussed 
further in Appendix B. 
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8. Heuristic arguments for the (non-) occurrence of isospectrality 

In the previous sections we have studied isospectrality for the cases s — \,N — 3,4 
and excluded certain other cases, e. g. complete isospectrality for s > \ and covariant 
isospectrality for N > 3. However, we have been unable to present a complete list of cri- 
teria for the (non-)occurrence of isospectrality What is also missing is some simple and 
intuitive argument why isospectrality is so rare. As a substitute for a complete theory 
we will, in this section, provide some heuristic arguments for the (non-)occurrence of 
isospectrality which also may give more insight into isospectrality than detailed proofs. 
We think that these arguments could be made rigorous as far as necessary conditions for 
isospectrality in the cases s = |,iV = 3,4,5 are involved. (See the remarks in section 2 
on the dimension of isospectral classes and the fibration theorem.) However, a detailed 
proof would require technical issues from the theory of algebraic varieties which are 
beyond the scope of this article and, moreover, would appear as superfluous given that 
isospectrality in some of these cases has already been proven by case inspection. 

The heuristic argument goes as follows: Isospectrality will (only) occur if the num- 
ber (^) of bonds between spin pairs exceeds the number L of independent eigenvalues 
of Hq. In this case, typically the systems corresponding to an n = (^) — L-dimensional 
sub-variety of J will possess the same eigenvalues. The argument may even be applied 
if one is not aware of all relations among the eigenvalues which determine L. In such a 
case of unknown relations one would perhaps over-estimate L and hence under-estimate 
the dimension n of the isospectrality classes but, depending on the case, one could cor- 
rectly predict the occurrence of isospectrality. 

When counting the number L of independent eigenvalues one first has to consider 
the (2S+ l)-fold degeneracy dictated by the rotational invariance of Hq. In the simplest 
example, one has to couple N = 2 spins s — |, obtaining one triplet and one singlet as 
eigenspaces of H , symbolically 2x2 = 3 + 1. Thus there are not four, but only two 
independent eigenvalues of Hq. Similarly, for iV = 3 one has 2x2x2 = 4 + 2 + 2, 
hence 3 independent eigenvalues. In the latter case there are = 3 bonds. Thus the 
heuristic argument does not yet explain isospectrality with n = 1-dimensional classes. 
However, it is easy to find a "missing relation" among the eigenvalues which reduces 
the number of independent eigenvalues to L = 2: It is just the relation Tr H Q = which 
yields a linear relation of the form AE 1 + 2E 2 + 2E 3 = between the eigenvalues E u . 

For s = | and arbitrary iV it can be shown that there are exactly (n^l) ^de- 
pendent eigenvalues due to rotational degeneracy and hence, considering Tr H = 0, 
L < (|_^y 2 j) — 1- I n this manner we obtain the results summarized in table 1. 

For large iV the the entries in the second column of table 1 grow asymptotically 
as 2 N J hence almost exponential, whereas (^) grows only quadratically. Therefore 
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Table 1. Occurrence of isospectrality for s = | and arbitrary TV based on a 
heuristic argument on the number of bonds (3rd column) and the maximal number of 
independent eigenvalues (2nd column). 



N 


([JV/2j) 1 


(?) 


Isospectrality expected 


2 


1 


1 


no 


3 


2 


3 


yes 


4 


5 


6 


yes 


5 


9 


10 


yes 


6 


19 


15 


no 


7 


34 


21 


no 








no 



Table 2. Non-occurrence of isospectrality for s — 1 and small N based on a 
heuristic argument on the number of bonds (3rd column) and the maximal number of 
independent eigenvalues (2nd column). 



N 


L < 


(?) 


Isospectrality expected 


2 


3 


1 


no 


3 


8 


3 


no 


4 


24 


6 


no 








no 



our heuristic argument will only predict isospectrality in the cases s = = 3,4,5 
but not for larger N. For larger s > \ the growth of the second column will prevail 
from the outset, see table 2, and our argument will not even be applicable for small N. 
Of course, this does not strictly exclude isospectrality for those cases, but makes it very 
unlikely in our opinion. 

There is another aspect which shows up in table 1 and which we now discuss: Note 
that for N = 4 the difference between the numbers in the 3rd and the 2nd column, 
6 — 5 = 1, would only explain one-dimensional isospectrality classes, whereas we have 
encountered two-dimensional classes in section 7. This indicates that the correct L 
should be 4, not 5, and that there is a further "missing relation" between the eigen- 
values of H , comparable to Tr H = 0. Indeed, as shown in the following paragraph, 
there holds a general property of the eigenvalues of arbitrary Heisenberg Hamiltonians 
H corresponding to the distribution of the eigenvalues among the quantum numbers S. 

Let, as above, denote by P5 the projector onto the eigenspace of S 2 with the 
eigenvalue S(S + 1). Then 

Tr(H F s ) = J*y ^(S x ■ S y F s ). (72) 

x<y 
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Since S 2 and hence all P5 commute with arbitrary permutations of spin sites, the last 
factor in ( [71]) does not depend on x, y and can be factored out: 

Tr (H F S ) = ( J2 J *y ) Tr ' »s) (73) 

\x<y / 

= JTr(S 1 -S 2 F s ). (74) 

Being independent of the J xy this factor can be calculated for any suitable H , e. g. the 
one with constant J xy = 1, 

H = ±(S 2 -Ns(s + l)), (75) 

which yields, after some computation, 

Tr(H F s ) = J j^—^s (S(S + 1) - Ns(s + 1)) TrP 5 . (76) 

Hence for all Hq, the vectors with the components (Tr(HoFs)) s=s s are 
proportional to the constant vector given by the r. h. s. of ([75|) with ,1 — 1. This 
gives a number of S max — S m i n = [Ns\ independent linear equations for the eigenvalues 
of Hq. Of course, Tr H = Tr(HoFs) = is a consequence of these equations and 
must not be counted seperately. For s — |, N = 4, 5 we obtain [Ns\ = 2 independent 
equations, which explains the two-dimensional classes for N = 4 we found in section 7, 
and predicts, at least, two-dimensional classes for the case N = 5 not yet analyzed in 
detail. If this case would show isospectrality classes of dimension n > 2 one could try to 
explain this by invoking more complicated relations derived from the higher moments, 
Tt{Hq¥s), k > 1. However, we will not further pursue this question here. 



9. Conclusion 



We summarize our results in table 3. It is in order to add some remarks on the 
possibility of determining the coupling constants J. Our results on the limits of uniquely 
determining the values of J in the case of isospectrality does not mean that these 
values could not be determined otherwise. First, we did not consider thermodynamical 
functions which do not come solely from the partition function, such as correlation 
functions, etc. Second, we do not adhere to a positivistic attitude which would in 
principle deny the physical reality of unmeasurable quantities. As in other domains 
of physics, these parameters could also be determined with the aid of additional 
assumptions, e. g. based on the symmetry of the molecules and supported by chemical 
considerations, which although plausible have not been confirmed directly. So we think 
the situation is different from theories with gauge freedom. 
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Table 3. Our results for the occurence of various types of isospectrality for different 



N i The 

1 V , O . -1 11C 


piUVLll ll^aLllLS all 


HHlll^ctljCd do JCS UI llu 




AT n 

IN S 


Plain IsospcctrE 


ility Complete Isospectrality 


Covariant Isospectrality 


3 ± 

2 


yes 


yes 


yes 


3 >i 


unlikely 


no 


unlikely 


4 i 

* 2 


yes 


yes 


no 


4 > i 

* ^2 


unlikely 


no 


no 


2 


likely 


likely 


no 


> 5 ± 


unlikely 


unlikely 


no 


>5 >i 


unlikely 


no 


no 
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Appendix A. Trace formulae 

Expanding the terms occuring in the trace 

Tr fl- = Tr J xy ]T Sg>S<A (A.l) 

\x<y i=l / 

and using Tr(A <g> B) = (Tr -A)(Tr B) one ends up with terms of the form 

Tr(Ai . . . At), A v e {S™, S^, S®}. (A.2) 

Here the spin operators without a site index denote operators operating in the single 
site Hilbert space C 2s+1 , not total spin operators. 

Let £i(i = 1,2, 3) denote the number of occurences of in the product A\ . . . At. 
One can easily show that Tr(Ai . . . At) is non-zero only if all £i are even or all £{ are 
odd. We give a list of the simplest cases, where the trace is non-zero: 

£ = : Tr(l) = 2s + 1, (A.3) 
£ = 2 : TiS^ 2 = ^ m 2 = \s{s + l)(2s + 1), (A.4) 

m=—s 

£ = 3: Tr S^S^S® = - Tr S^S^S™ = l -s(s + l)(2s + 1). (A.5) 
From this we obtain for n = 2: 

x<y i=l 
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= (E^)r 2 ( s + 1 ) 2 ( 2s + 1 ) JV - (A-7) 

\x<y / 

For n = 3 there occur two kinds of non-zero terms: 

Tr {S^S^S^ ® S^S^S® ®l w _ 2 ) (A.8) 
and those terms obtained by permutations of {1, 2, 3}, and 

Tr (Si t)2 ® S« 2 <g) S« 2 ® Iat-s) , (A.9) 

where x, y, z are pairwise distinct and % = 1,2, 3. Consequently we obtain: 

Tr Hi = T X + T 2 , (A. 10) 

7\ = ^43!(2s + l)^ 2 (Tr 1 S« 1 S( 2 )^ 3 )) 2 (A.ll) 

= -E J ^ s2 ( s + 1 ) 2 ( 2s + 1 ) iV ' ( A - 12 ) 



'6 

:r<{/ 

T 2 = 3! ( J2 J*yJ* Z JyA (2s + l)^ 3 ^ (Tr 5« 2 ) 3 (A.13) 

\x<y<z / i=l 

= ( J2 J *y J *z J yA \s\s + lf(2s + l) N . (A.14) 

\x<y<z / 

Appendix B. The exceptional cases 

Here we collect some properties of isospectral classes for N = 4, s = 1/2 which belong 
to the exceptional case. Although we did not obtain a complete classification of all 
isospectral classes these results may be useful for further studies. 

The subset J\TZ, which denotes the complement of TZ, is characterized by the 
vanishing of all 4 x 4 sub matrices of hence it will be also an algebraic variety. After 
some computations one shows that J\1Z is a union J\1Z = C\ U C 2 U C 2 of three simpler 
varieties given by the following equations 

d : h = D\ 2 + D 2 13 + D 2 14 = 0, (B.l) 
C 2 :I 2 = 3/ 3 2 , i. e. S 12 = S 13 = S u , (B.2) 
C 3 : / = 9 = 0, (B.3) 



where 



/ = {D\ 2 - D 2 14 )D 13 + J D 12J D 14 (.S 1 4 - S 12 ) = 0, (B.4) 
g ee (Dj 3 - D 2 U )D 12 + D 13 D U (S U - S 13 ) = 0. (B.5) 



In the first case it is clear that an isospectral class lies entirely inside C\ or outside 
Ci, since ii is constant on this class. An analogous remark applies for C 2 but the case C 3 is 
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more complicated. In the first case with /i = the isospectral classes are O-dimensional, 
since the circles defined by I2 = H and ^3 = i% intersect the variety I4 = —S12S13S14 = 14 
at most at six points. In the second case with S12 = S13 = S14 = a the isospectral 
classes are the 1-dimensional intersections of the sphere I± — > and the variety 
I4 = 2Di 2 Di 3 Di4 + (Df 2 + Df 3 + Df^a — a 3 = 14. In the third case we have shown, using 
the Eliminate-command of MATHEMATICA® , that / = g = implies an equation of 
the form P(I\, h, h, h) = 0, namely 



Hence those isospectral classes with P(Ii, I2, I a) ^ 0,/i ^ 0,^2 7^ 3/f lie entirely 
inside 1Z and belong to the generic case. We conjecture that also the converse holds, 
namely that P = implies f — g — but could not prove it. By calculating the 
corresponding Groebner bases, it can be shown that P and /, g generate different ideals 
in the ring of polynomials in 6 variables. However, this will not exclude the possibility 
that the corresponding real algebraic varieties may be equal. Hence we cannot exclude 
further exceptional cases in the realm P = but / 7^ or g 7^ 0. The isospectral 
classes in the case / = g = will be 1-dimensional: The constraint / = g — allows 
one to express the variables S12, S13, S14 in terms of I 3 , D 12 , Di 3 , D u . The remaining 
constraints h — ii, I4 — 14 define a family of curves obtained as intersections between 
spheres and tube-like surfaces. 
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